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SIHIARY 


This paper discusses the NASTRAN experiences of a major structural 
design and fabrication s\d>contractor that has less engineering personnel 
and computer facilities than those available to large prime contractors. 
Efforts to obtain sufficient coB|>uter opacity and the development and 
implementation of auxiliary programs to reduce manpower requirements are 
described. Applications of the NASTRAN program for training users, diecking 
out auxiliary prograats, performing in-house research and development, and 
structurally analyzing an Avco designed and manufactured missile case are 
presented. 


INIRDDUCTION 


The Avco Aerostructures Division has long been actively engaged in the 
researdr and developamnt areas of structural analyses. Since 1966, particu- 
lar emphasis has bein placed on finite element techniques utilizing the 
dinlacement method* Sudi efforts yielded a static analysis program for 
bow the IBM 360/Model 40 and IBM 1130 conputers* Although the program was 
someidut limited in the types of finite elements, the objective for solving 
structural problems containing a large nunber of degrees of freedom was 
achieved. This objective, however, was not obtained without the sacrifice 
of coBputer efficiency due to a requirement for large anoints of peripheral 
processing time* For eocaqple, the IBM 1130 could acoonaodate a structural 
model with nearly 1000 degrees of freedom, but approximately 10 hours were 
necessary for a solution* Run tines of 3 and 4 hours were not incommon for 
IBM 360/40 analyses* Only bar, triangular meahrane plate, and rectangular 
amabraM plate elaMnts were available in the IBM 1130 prograa, but the IBTf 
360/40 prograa also included the triargular bending plate element. (Plate 
elements assumed constant stress conditions*) It becaw immediately qiparent 
that the coagilete cipabili^ for performing structural analyses (i.e. , static, 
vibration, buckling, etc.) by finite element tedniques was not within the 
practical i^iper limits of Avco A/D's coa^uter facilities* 

In the latter part of 1970, Avco A/D became aware of the Industry 
Researdi Associate Prograa initiated by the NASA-Langley Researdi Center* 

The value of the program for providing a mutual interchange of technology 
was iaaediately rao^nized and Avco became an active participant in January 
1971* The initial participation consisted of two engineering personnel 
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being assigned to the Structitres Division of the Stmcture Directorate* It 
was in this time period that the NASTRAN Systems Management Office (NSMD) was 
established in the Structures Division at NASA-Langley (October 4, 1970) and 
the first public release of NASTRAN through OOSMIC occurred (November 1970) . 
Since NASTRAN was being heavily used at N^A-Langley and other government 
centers for testing and evaluation, NASA- assigned Avco personnel could evalu- 
ate the program and beccme qualified users v4iile performing their assigned 
tasks* This NASTRAN ei^rience is one excellent example of the spin-off 
benefits from the NASA/Industry Research Associate Program. 

It was obvious from the initial confrontation of the Avco associates 
with the NASTRAN program that it would enable us to attain our research and 
development goals. However, the in-house cooqauter capacity at Avco A/D \;as 
still not adequate to accommodate the program* A survey of the local 
(Nashville) conq;)uting facilities was undertake to determine if any conputers, 
suitable for NASTRAN operations, were available* In-house terminals, as 
supplied by conputer leasing services, were disregarded since the NASTRAN 
usage was upredictable; operating costs could be excessive on a per run 
basis. A colter service, NLT Computer Services Corporation, which primarily 
serves the local banking and insurance interests, was fotnd to have centers 
conpatible with the NASTRAN program; namely, an IBM 360, Model 65 and an IBM 
370, Model ISS oonputer* Arrangements were made with NLT for implementation 
of the level 12 NASTRAN on each of their conputers* Further, it was decided 
that, vdienever possible, auxiliary NASTRAN programs capable of generating and 
checking input data and manipulating ou^ut data (i.e* , resizing routines) 
would be incorporated on the in-house IBM 360/40. The inconvenience of the 
data management problem presented by this approach was considered less 
indesiral)le than tlie economic problem of using the NLT conputer services for 
all cases* Avco A/D acquired the NASTRAN program in October 1971 and had it 
operational on the NLT conputers qiproximately one month later* The inple- 
mentation of the program by NLT personnel will be discussed in this report* 
Further, this unipie arrangement between the financial and technical connunity 
will be described with particular enphasis on the problems encoxmtered due to 
different terminologies and concepts. 

A description of Avoo's auxiliary conputer routines to be used with the 
NASTRAN analysis program is included in the paper* For exanple, the data 
generation programs described in references 1 and 2 have been converted from 
a CDC 6000 series conputer to the IBM 360/40 conputer with some modifications* 
Also, a program vhidi checks the NASTRAN iiput data for format and syntax 
errors and inconplete and/or duplicate data and generates a tqw for plotting 
undeformed structure lias been di^loped* This dwck-out program is s&iilar 
to the special NASTRAN program described in reference 3* 

In addition, results obtained from various NASTRAN investigati<ms of 
widely varying t]^s of structure are presented* For each of the investi- 
gations, the results describe (perational problems, run times, and core 
requirements (including oonparisons of the IBM 360/65 and IBM 370/155 com- 
puters whenever possible) plus conparisons of the NASTRAN results %dth 
theoretical or eiqperinental data* 
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The pq>er also discusses the in^rtance of a NASTRAN type program to 
Avco Aerostructures Division as a major subcontractor to the prime aerospace 
contractors* 


AVCO/NLT NASTRAN IMPLEMENTATION AND OPERATION 


An interfac^ of the capid>ilities and facilities of the Avco Aero- 
structures Division and the NLT Gonputer Services Corporation was required 
to acccmnodate the NASTRAN program. The acquisition (and ^dating) and 
execution of NASTRAN was the responsibility of Avco* The inplementation and 
maintenance of NASTRAN was the responsibility of NLT* Obviously, NASTRAN 
bridges the gap that exists between the ci^abilities of the technical and 
financial commnities; namely, scientific progranners* 

Before NASTRAN, NLT had provided a data processing service only and was 
con|>letely ine)q;>«rienced with respect to scientific con;)uter programs. How- 
ever, the program was ijiplemented in approximately one week and operational 
in less than one month on the IBM 370, Model 155 and the IBM 360, Model 65 
conputers* The primary difficulty encointered in the in^lementation of the 
program was NLT^rsonnel's lack of familarity vrith the OS operating system* 
(NLT uses the DCS operating system for their data processing services*) Minor 
problems occurred ^ to a reluctance to believe the disk space requirements 
of NASTRAN* Some check-out problems were encouitered in the execution of 
the demonstration problems because of core size requirements* Further, the 
pinched output for the restart deck of demonstration problem 1-1 was com- 
pletely unexpected* Q[>u times for the execution of certain demonstration 
problems on the IBM 370/155 conputer are presented in table 1 along with compat 
ible cpu times obtained from various computers as given in reference 4. 

Although more specific and concise docunentation pertaining to the inplemen- 
tation would be desired, the ease in making NASTRAN operational has been 
substantiated* In all cases, COSMIC was very prompt in identifying user prob- 
lems and giving a solution. 

Some operational problems have been encountered due to the conflict of 
NASTRAN with the previous experience of NLT personnel * A limit of 900 cpu 
seconds was a standard NLT cooputer exit* This limit was based on previous 
eimrience vhich normally indicated the presence of a progranming error (i*e*, 
a 'hard DO loop') • Execution of NASTRAN has made them accustomed to exceeding 
the 900 qiu seocnds restriction vMch is now deleted front all NASTRAN prob- 
lems* The computer operator confused the OFTP (old problem tape) and NPTP 
(new problem tape) with the outyut and iiput tq^s, respectively* Whether 
or not this contributed to an apparent pimlem of tqpe management was never 
identified. 

It is believed that this arrangement permits smaller engineering depart- 
ments to attain a structural analysis opuility that may otherwise be 
iapractical and uneconoadcal. 
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AUXILIARY PROGRAKS 


Since the different physical locations of the Avco engineers and the 
NLT conputers presented a data handling problem, it was decided that NASTRAN 
auxiliary programs would be installed on Avco A/D's in-house IKI 360, Model 
40 conputer. Further, those routines which would check hand generated input 
data or automatically generate the input data offer an economic advantage in 
the decrease of aborted analyses (cominiter cost) due to input errors and the 
reduction of manpower for data preparation. 

A special NASTRAN program which checks the input data for mispunched 
data cards and incorrectly transcribed data and goierates a tape for unde- 
formed structure plots has been developed for the IBM 360/40 corputer. This 
program is functionally equivalent to the check-out program for the OX! 6000 
series omputer described in reference 3. The program operates under the 
05 360/PCP operating system and requires 210K bytes of a 256K byte core and 
a disk drive for perijmral storage. Only ‘ sii»routine BTSTRP required a 
revision due to the necessary reduction in the length of a GINO buffer from 
1803 to 250. This change caused a corresponding reduction in the default 
parameters of the block sizes contained on the PROC job control card to 
SP2-1, BLKl-1028, and BUC2-1032. 

The automated input data generation routines, as described in references 
1 and 2, have been converted from the QIC 6000 series conputer to the IBM 
360/40 computer. These routines offer a siibstantial reduction in the man- 
hours required for model generation since they recpiire a mininun of input 
^ta tdiidi can be obtained from engineering drawings. This manpower reduction 
is particularly attractive to smaller engineering departments such as Avco A/D 
because a cotplex structural analysis can be perfoxm^ by a minimun nunber of 
engineering personnel* In conducting some recent NASTRAN analyses, finite 
element models of a delta wing and a sequent of a missile case were generated 
by these automated routines in less than 4 man-hours each. The generality of 
the routines was demonstrated in reference 2 and their versatility has per- 
mitted the modeling of a railroad passenger cat by Avco A/D. Another routine 
vhich automatically generates a triangtslar or qua^ilateral mesh about a 
circular cutout in a flat panel, as investigate in reference 5, was modified 
to yield punched NASTRAN input «te. Presently, programs are being developed 
which will calculate the stiffiiess coefficients for ^cialized structiural 
eleeants (i.e., reinforced concrete or inUgrally stiffened plate elements) 
and pmdi the data on cards wi^ an input foxmat consistent with NASTRAN. 

The contour plotting routines that are described and demonstrated in 
reference 1 have been converted from the Q)C 6000 series coaputers to the 
IBM 360/40 coaputer. Since Avco A/D lacks an in-house plotting capability, 
the plotting routine has been modified to yield a printed output of the 
planar coordinates of each contour line to peimit manual plotting. A plotting 
service, similar to the Avoo/NLT coaputer amngenent, could be obtains since 
plotters that are coapatible with the routine md NASTRAN are available in the 
Neshville area. 
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Currently, a fully stressed design techniaue is being developed vdiich 
is almost identical to the routine incorporated in the " program (ref. 6) 
and the SAVES program (ref. 2). The input and output t;.iT.'>at' ^or this routine 
are c^atible with the NASTRAN output and input for ats, res^crtively. A 
simplified design criterion for the sizing of rod and bar elements has been 
incorporated into the program. 


APPLICATIONS 


Since the acquisition of NASTRAN by Avco A/D in October 1971, it has 
been used for perfotroing in*house researdi and development programs, check* 
out and verification of auxiliary routines, and analysis of existing 
structural designs* The following discussion is directed to three investi- 
gations of widely varying types of structure* 

Reinforced Concrete Slab - The purpose of this study was to evaluate a 
method of modeling reinforced concrete with finite elements that are suffi- 
ciently accurate and economically feasible. Quadrilateral membrane and 
bending plate elements were used to represent the concrete slab and rod 
elements were used to model the reinfoi'cing rods. The representation of each 
reinforcing rod *#ould yield a finite element model with a maximun mesh size 
dictated by the reinforcement spacing. Since this pacing is relatively 
small, it is conventional to represent several reinforcement rods by a single 
rod eleaient which is positioned between two nodes along the edge of a larger 
plate element* This method is coanonly called the lunping technique and 
results in large rod elements at the neutral axis of m plate's cross 
sectional area. However, the reinforcanoit rods are usually displaced to the 
tension side of the neutral axis and thus cannot be included by the conven- 
tional luMping technique. For this study, continuity of the reinforcing rods 
and concrete (no sliming) was assuaed. The elongations of the rods are 
considered equal to mir noimal displaoeamnts on the edgewise face of the 
plate element at the intersecting rods* A linear variatim of the trans- 
lational md rotational displacements is assuwd between two nodes of a 
particular edge for conputiM the normal displaceaients* This coagnitation of 
the normal displaceamnt at the reinforcement rod position was accoeplished 
by use of the multipoint constraint equations contained in NASTRAN. 

The subject of the irivestigation was a square, simply supported, rein- 
forced concrete slab with reinforcing rods spaced at 2S.4 cm. intervus along 
each 609*6 cm. edge. The slab was 20.32 os. thidc md si;bjected to a dead 
weight loading. The roinforemnent rods (diameter • 1*90 cm.) were at a depth 
of 3.49 cau from the tension side of the concrete. The non-linear properties 
of the concrete due to its inability to withsund tensile stresses were 
neglected. Model mesh sites of 4 x 4 md 6x6 (corresponding to element 
sites of 152.4 X 152*4 SS; mnd 101.6 x 101.6 cm*, respectively) were investi- 
gated* Figure 1 presents a ooo|)arison of tOie rod stresses for eadr NASTRAN 
malysis with the theoretical results determined from closed form equations 
on page Vi6 of reference 7. Figures 2 mid 3 show a coeparison of the maximua 
comiressi'O stresses in the concrete mnd lateral displacement of the slab, 
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respectively* Whenever necessary* NASTRAN results are extrapolated to the 
desired locations* 

Ihe coarse mesh model (4 x 4) required 218 cpu seconds for a solution 
on the IBM 370/155 computer and 260 cpu seconds vere required for an IBM 
360/65 computer solution of the fine mesh model (6 x 6) using Level 12 
NASTRAN* Due to the difference in problem site* it was not possible to 
evaluate the confuting speed of the different centers* However* the 
prohibitive time required for processing and manipulating multipoint con- 
straint relationships was indicated by this investigation, i^roxiniately 
201 of the total q>u time was used on the multipoint constraints for each 
model* It is because of this time penalty that the previously mentioned 
auxiliary program for generating NASTRAN pinched input of stiffiness coef- 
ficients for specialized structural elements is being developed. The 
NASTRAN run time will be further reduced by the time required to calculate 
the stiffimess* Obviously* the NASTRAN output would have to be used in an 
additimal auxiliary program to determine the actual internal loads and 
stresses (i.e** axial stress in reinforcement rods)* 

Ihese preliminary NASTRAN results show good agreement with theoretical 
results* Further validation of the computational efficiency and accuracy is 
required when the auxiliary routines become operational. This metliod is 
also applicable to integrally stiffened or semi-monoooque panels and other 
structures vdieiein an assemblage of finite elements is represented by a 
discrete structural element* 

Structural Effects, of W ing Vfrittiflng ‘ 'Hiis study was primarily 

initiated to check out the lifting surface data generation routine* To 
provide additional NASTRAN experience* it was decided to conduct a siirple 
design stady to evaluate the change in stiffiitess and strength due to camber 
variations of a wing* The %ring has a delta planform and a modified diamond 
airfoil section wiu a 3% maximuR thiclousss* A uniform pressure loading was 
mlied to the structural ludel* For a maxi min camber of 1*501 chord* the 
finite element model is identical to the luqped finite element model of the 
baseline wing structure used in the study described in reference 2* The 
normalized tip deflection and stress distribution in the lower surface skins 
near the wing-fuselage intersection fbr varying amoimts of camber are shown 
in figure 4* 

The finite element amdel of each %dng configuration is identical except 
for the 'out of pUne' geometry dmnges caused w camber variaticais* If 
everything r ama ina the saam eacoapt for wing oonfiguretion* the qpu time 
should be constant* Althomh there %«s soma tmeiqplained but minor variations 
in the cpu times* approadmately 490 and sao aeconds were required for NASTRAN 
Level 12 solutions on the IBM 360/65 and IBM 370/lSS computers, respectively* 
These times are contradictory and have not been e)q>lained by NLT systems 

K rsonnel except to speculate that it could be due to the operating system, 
all cases* the 06/MVT/HASP operating system was used* 

JilMXSRiri? * This stu^ was initiated to detensine a possible cause 
for observable damage on Avco designed and manufactured missile case sections 
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due to operational trsisportation conditions, llie objective of the investi- 
gation was to determine the cause of the damage and study possible design 
changes in the transporter to prevent the re-occurrence of structural damage. 
Ihe sections of concern are lightweight honeycomb shells with ring frames at 
the ends of each section. The damage occurred at one of the transporter's 
siq^rt saddles which was located such thct the honeycomb shell would rest 
on it. The missile case sections were modeled by OQUADI elements for the 
shells and OAR elements for the frames. The loads apx>lied to the models 
were in accordance with loading conditions used for testing* Clanged boundary 
condi'wions were assumed at the intersection of the honeycomb structure and 
the relatively thick-walled back-v|) structure. Those grid points in contact 
with the saddle were not allowed to translate in the lateral direction. 

Results from a correspond^ analysis and the test condition are shown in 
figure S in tdiich the variation of the outer face sheet longitudinal (membrane) 
stress about the circumference of the shell is presented. The variation of 
the test data (* e and - e) is due to the lack ox symnetry of the structure 
caused 1^ cutouts for doors and access ports. Although a good correlation of 
NAST^ and test data was obtained, some variations should be expected since 
the cutouts were ignored in the finite element model* However, a NASTRAN 
analysis did predict a core shear failure at the exact location v4\ere 
observable danuige did occur as indicated in figure 6. 


REOMCNnAriONS AND OONOUSIONS 


Based on the NASTRAN ejqMrience obtained to date at Avco A/D, the 
following chflfiges or improvements would be desired. 

1. The development of a mathematical technique for estimating the time 
and core requirements for finite element <*H)dels of various sizes md 
degrees of coeplexity would reduce turn aroind tine. This technique 
could be in the form of ee|)irical relationships amenable to pro- 
gnmming on an office size ooe|>uter. 

2. Since the static analysis portion of NASTRAN receives the greater 
usage, a mini NASTRAN which only contains the rigid format 1 and 
reduced buffer sizes to allow for execution ui srialler computers 
(i.e*, IBM 360/40) would probably prove economical. 

3. The margin of safety pertaining to ooBf)ression members i . not con- 
ventional ainoa a positive margin of safety is under strength and 
vice versa for a NASTRAN analysis* This msy be corrected by using 
a negative m»ber xtt designate the allowable comf ressive stress ot> 
the materials p - aperty card* 

In conclusion, the end continuing maintenance and improve- 

ment of the NASTRAN pregrstt has made it possible for coipanies with limited 
engineering manpower and ocaputer facilities to attain a proficiency in 
at uctural analysis. Prior to NASTRAN, a large amount of the independent 
research and de^lopam nt ftsids was expe n ded on the developner>t of structural 
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analysis ..lethods based on finite element techniques. This effort can now be 
directed to structural research utilizing NASTRAN as the analysis tool, 
^rther, the inplementation of NASTRAN by systems personnel v^o are completely 
inexperienced with scientific programs and comprehension of the program by 
new userw attests to its usability. 
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TABLE 1. OTITARISON OF LE\T;L 12 EXEaTTIGN T»ES (0>U SECDNES) 
FOR NASTRAN DRDNSTRATION PROBLEMS 


Gon^uter Ccnfiguraticn/Series/Model 


PTOblCTI 

Fom § 

IBM 

IBM 

IBM 

CDC 

UNIVAC 

No. 

Version 

370/155 

360/95 

360/67 

6600 

1108 

1-1 

U 

116 

60 

108 

31 

72 

I-IA 

R 

72 

30 

102 

24 

24 

1-2 

U 

166 

84 

264 

72 

60 

1-3 

II 

242 

114 

336 

120 

96 

1-4 

U-1 

1022 

020 

1464 

666 

1020 

1-4 

U-2 

2665 

380 

-- 

■» A 

• « 

1-5 

II 

1505 

660 

1800 

m m 

1080 

1-6 

U 

68 

48 

132 

30 

24 

1-7 

U 

251 

60 

270 

120 

90 

2-1 

U 

88 

72 

168 

30 

24 

3-1 

U-1 

1567 

720 

3900 

826 

660 

3-1 

U-2 

3088 

610 



• w 

4-1 

U 

328 

168 

486 

144 

90 

5-1 

U 

863 

408 

1134 

486 

318 
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NORMALIZED AXIAL STRESS 


NOTES 1. STRESSES AT Ya 0, NORMALIZED ON 
MAXIMUM THEORETICAL AT X»0 
2. ROD SPACING* 25.4 CM, 

SEMI-SPAN, A- 304.8 CM. 


1.00 


0.80 


0.60 



Y 




0.40 


0.20 


0 



□ THEORETICAL (REF. 7) 
O NASTRAN 4X4 GRID 
V NASTRAN 6X6 GRID 


20 40 60 

PER CENT SEMI-SPAN 



Figure 1. • Variation of rod axial stresses in a square i 
singly supported* reinforced concrete slab 
si^jected to a dead weight loading. 


604 



NORMALIZED COMPRESSIVE STRESS 






NOTES: 1. STRESSES AT Y»0. NORMALIZED ON 

MAXIMUM THEORETICAL AT X»0 
2. SLAB THICKNESS-20.32 CM. 
SEMI-SPAN, A-304.8 CM. 

Y 



Figure 2. * Verietion of concrete oonpresslve stresses in a square 
sisply siimorted, reinforced concrete slab 
subjected to a dead v«ight loading. 


NORMALIZED LATERAL DEFLECTION 


NOTES: 1. DEFLECTION AT YsO, NORMALIZED ON 
MAXIMUM THEORETICAL AT X-0 
2. SEMI-SPAN, A»304.8 CM. 



Figure 3« - Variation in lateral deflection of a square, 
siaiply siyiported, reinforced concrete slab 
subjected to a dead weight loading. 
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Figurt 4» • Strtfs dictributian near wing-fuaelage intersection 
and nonalized tip deflection for vatyutg degrees of canber. 
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LONGITUDINAL STRESS, KN/M 


o TESTDATA+e 

□ TESTDATA-0 
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Figurt 5. - Variation of outer face sheet longitudinal stress 
at M.S* 70. S for synmetric vertical loading. 
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Figurt 6. * CiroiRferential transverse shear flow (core) 
for cradle sv|>port enclosed angle • 90* 
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